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EXAMINATION OF THE VALIDITY OF.AN APPEOXIMATE 
SOLUTION OF A CEETAIN VELOCITY EQUATION. 


By a. Brown, M.A., B.Sc. 
(Eead November 28, 1907.) 


§1. This problem has been suggested by the important part 
played by “ initial conditions ” in the solution of differential equa¬ 
tions occurring in Dynamics. In certain cases the special character 
of the “ initial conditions ” makes a simple treatment of the equa¬ 
tion possible ; but sufficient attention has not been paid to the 
question whether the simplification introduced is legitimate, or to 
the investigation of the limits within which it holds. 

To make the matter clearer, consider the equation— 


dy 

dt 


-y+ky^, 


with initial condition y = \ when i = 0, X being supposed small. 

One method of obtaining an approximate solution is to start with- 


dy 

di 


= -y 


and find tbe solution of this equation. The result is 

It is then assumed that if the solution thus obtained is small for 
all values of t, the result given is an approximate solution of the 
original equation. 


The actual solution of the equation is X . e' 


and in 


• 1 - XA:(1 - £-‘)’ 

order that this be expansible in powers of X for all values of t it is 
necessary that Xk<l. If XZ!:>1 the solution of the simplified equa¬ 
tion is not even approximately a solution of the original equation ; 
and unless Xk is a small fraction the solution given will not be a 
good approximation. The validity of the simplification depends on 
the nature of the coefficients as well as on the smallness of X ; and 
it seems important to investigate how small X must be for the 
simplified form to be used, 
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§2. The first question we consider is— 

Under what conditions can a solution of— 

. ■ .( 1 ) 

with initial condition y = X when i=0, be obtained in a series of 
ascending powers of X ? 

[This simple form is taken in order to illustrate the method used 
in the more general case, and also because the result can be tested 
by comparison with the solution in finite terms.] 

To determine the possibility assume a solution of the form— 

y = T„+XT, + X^T,+ .(2) 

and determine the T’s so as to satisfy the differential equation. If 
the series thus obtained is convergent the existence of a solution of 
this form is established. 

From the initial conditions we have— 


To = 0, T, = 1, T 2 = 0, T 3 = 0, ... when ^ = 0. 

Substitute from (2) in (1)— 

.•. t;+xt;+x-t:+...= a.(T,+ xt,+x^t^+...) 

+^^2(T^^-XT, + X-T2...)^ 
dT'’ 

where T[. means . 

Equating powers of X— 

T; = a,T,+a2.2ToTj I 

T; = a,T2+a,(T;+T,.T2)| . 

&c. > 

The first of these has the initial condition To=0 when ^ = 0 ; hence 
initially T[=0 ; and thus To=0 for all time. 

Transform the set (3) by the substitution— 

T.^S^.e^^xt.(4) 

and the equations (3) become a set of which a typical member is— 

s;=a2. £nSoS„ + S,S,_.+ ... + S,_,. S, + S,So] .(5) 

The initial condition applicable to Tj gives us Si = l when ^=0. 


Further substitute— 


s,,=ar'E, 


(6) 
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for values of r above unity and the typical number (4) becomes— 


E,.—[EoEy+E iE2_, .. • + Rj_iRo] 
Hence, using the initial conditions— 


(7) 



It is hei'e to be noted that E,o=0; and hence when the expansion (2) 
is applicable the values of the coefficients can be determined by 
quadratures from (8). 

The solution is— 


(9) 


Xfrtii [1 + fijA-Rs + U2\^R3 + ...] 


To discuss the convergence of (9). 

From the values of Rq and Rj, and the mode of formation of the 
subsequent coefficients, it is clear that the R’s are all positive for 
positive values of t. 

is one factor in the solution, and will occur to higher powers 
in the other factor; so that, if we include in our consideration large 
values of t, a necessary condition for the convergence of the series is 

negative. 

Put Ui= — a, where Uj is positive. 

Applying this result we have— 


Ro=0 
R, = l 



rt 


^(RjRj + R2 Ri)< — (KjKj + K.1K2) 


( 10 ) 


10 



■n—i ( 


) 


where the K’s satisfy the relations— 

K, = 0 
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This can be simplified still further by putting— 


(excluding r=0. 


We now have— 


y. 


Lo=0 

L, = l 
L,=:l 

L'3 = (L3Li + LiL,) 

(^n—“I” • • “1~ —2 ”l~ i) 

To determine the L’s put— 

so that l+a;^f = L+L 2 a; 


ic2(l+= (Lo+L^a?+ Ij,x^ ...y 

= Lo+ar(LoLi4-LiLo)+;r^(LoL,+LjLi + LjLo) +.. 

=x^z 

{iA-x^y=^ 


_1—a?2— ^1—4a; 


, selecting the appropriate sign 


= ^c,+2C3.a7+2.c,a:^ + . 


where c,. is the coefficient of - x^ in the expansion of (1 - 4a;)^. 

Hence L,.=^c,.. 

Now the expansion of (1 —4a?)'' in ascending powers of a? is conver¬ 
gent if a? <\; hence and L^^-:^'4. 

^ n—cc 11=00 


11=00 


Thus the solution of the differential equation is— 


Xe 


where 4. 


.r., a^X _-r aA^-r-r "1 

_a.f .H,+ -^H3... 

L Oi a, ^ J 


This will therefore be convergent if 


uA 


<4- 
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Hence conditions which are certainly sufficient for the applicability 
of this form are — 

( 1 ) ai negative. 


(2) X< 


4a, 


(For the case considered in §1 this gives whereas we know 

that the range of availability of the method is bounded by 
§3. Apply the same method to the more general equation— 

^=a^y-\-ayj^-\-a^y^+ .(13) 


subject to the initial condition y = \ when ^ = 0. 

Try a solution of the form— 

2/ = To + \Ti + X2T, +. 

Substitute in (13) and equate coefficients of X— 

•'•To = aiTQ + a,T'+ a 2 T^ +... 

T'l=aiTi + a.2(ToTi + TiT^) + a^^ . dT^Ti + ... 

Tj = ttiT, + a,(ToT, + Tj) + a^. 3T^. + — 

&c. 


(14) 


(16) 


The initial conditions for the different T’s are— 

To=0, Ti = l, Tj = 0, T 3 = 0 ... when ^ = 0. 
Substitute Ty=£*i^S,.— 

So =0 

S. = l 

SU«x*. a,(SoS, + S?) + . a3(3SoS0 + ... 

S; = £«^«. a,( 2 SxS,) + £-«xt(3SoS3 + GSxS.S^) + ... 


(16) 


Once more it is to be noticed that the coefficients T or S can be 
immediately evaluated by quadratures, since the coefficient of T, in 
the equation for T^ is zero, and the earlier T’s are calculated in 
succession. 

As in the case discussed in § 2 , a first condition for convergence 
will be ai negative, equal to - ai (say) where aj is positive. Follow 
the line taken in § 2 , putting \ a. \ =a., |a 3 | = 03 , &c.; we then find 
that each of the S’s has, for all values of T, a modulus not greater 
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than that of the corresponding E where the set of E’s is given by 
the equations— 


Write- 


Eo =0 

E, = l 

E, = “^EoE, + El) + ^(3E„E1) + .. . 

cti 

E 3 = ~ . 2ExE, + ^ ( 3 E 0 E 3 + GE^EaEo) + • • • 
&c 

1 } = Ei^ + Ej^^ + E34^ + ... 


.(17) 


a- 

Then — . ip + Q — . + n . if + ... 

a, 2ai Octi 


= —^(Ei^ + Ej^' + .. .)^ + 2 J (Eli' + E,^^ +.. .)3 + ,.. 

= E, 4 ’ + Eoi;^ + + ... 

= 1] - Ei4 


• 'C 

• s 


= n 


a, 

— • V' 
«! 




The reversion of this series gives an expression for in ascending 
powers of 4 , the coefficients therefore being the E’s. 

In a previous paper * I have shown that the reversed series is 
convergent if U| <J, where J= ^(l + 2/3)" + l —(1 + 2/3), /3 being 

the greatest of the set — 4 ^... 

^ a/ 2af 3ar 

Hence finally the solution of —■ 

dy 

-^=a,y + anj- + a^y--]r... 


with initial condition y = \ lohen t = 0 can certainly he obtained by 
successive approximation in a convergent series of poioers of X pro¬ 
vided 

a-i is negative, 

and IX I < ^1 + 2/3)' + 1 — T ^/^)» 


'where I3 is the greatest of the set 




an 
_ ^ 


^4 

a. 

j 

2ai 

j 

3ax 
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